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Abstract 

An interplay of optical lattices and nonlinear impurities in controlling the dynamics of Bose- 
Einstein condensate bright solitons is investigated using effective potential approach. The ability of 
pushing the solitons into or away from the impurity region by changing both lattice and impurity 
parameters is suggested. A possibility for the existence of stable fundamental gap solitons, which 
appear to satisfy an inverted Vakhitov-Kolokolov criterion, is examined. 
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1 Introduction 

The first experimental observation of Bose-Einstein condensates (BECs) in 1995 [T] is one of the great 
achievements of modern science. This milestone achievement has been motivated many theoretical and 
experimental works j2j including creation and propagation of matter- wave solitons [3]. Solitons in a 
BEC are traditionally produced by changing the sign of atom-atom scattering length (a<j) from positive 
to negative values or vice versa with the help of magnetic |1] or optical [5] Feshbach resonances (FRs). 
These FRs either individually or a combination [6] of them offer an additional flexibility to obtain spatial 
variation of the scattering length and thereby introduces spatially inhomogeneous nonlinear interaction 
in the condensate. BECs with such nonlinear interactions have become systems of considerable current 
interest since it is expected that these systems may help to understand the dynamics of solitons emission 
which will be useful to construct atom laser [3 E]- In view of this, some theoretical works have been 
done on the existence, stability and localization of matter-wave solions in single -component O |TOl [11] 
as well as two-component [121 [131 E] BECs with spatially spatial inhomogeneous atomic interaction. 

A BEC can also be loaded in periodic potentials. These potentials, often called optical lattices 
(OLs), are generally created by interfering two counter propagating laser beams [I5]. The geometry and 
depth of the OL can be completely controlled [16] to create a large variety optical potentials. Studies 
in the dynamics of BEC solitons in periodic potentials with different spatially inhomogeneous atomic 
interaction have received a great deal of attention. Depending on the spatial variation of the atomic 
interactions, this system shows some new and interesting phenomena. For example, a localized state of 
BECs loaded in OLs can undergo delocalized transitions [17J when the atomic interaction is made to 
vary periodically in space. This inhomogeneous interaction helps a BEC in OLs to support gap solitons 
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in the upper half of the forbidden gap even when oq, the background atomic scattering, vanishes [18]. A 
special feature of a bound state in a BEC with localized inhomogeneous interaction is that it can contain 
any arbitrary number of solitons |H]. Moreover, the potential arising from the spatially inhomogeneous 
atomic interaction may help a BEC to support bright as well as dark soliton solutions [I9]. Recently, 
Sakaguchi et al investigated the existence, stability and mobility of solitons in single-component BECs 
for spatial periodically modulated nonlinearity [20]. For two-component BECs studies on the stability 
properties of solitons in OLs with spatially periodic interatomic and inter-species interactions have also 
been made in [2T]. The object of the present paper is to study the dynamics of matter-wave solitons 
in OLs with spatially localized (Gaussian) nonlinear interactions within the framework of variational 
approach. In particular, we work with a cigar-shaped BEC and make use of a potential model to 
investigate effects of the optical lattice wave number as well as width of the Gaussian nonlinearity 
on localization of the BEC soliton. In addition, we examine that the stable fundamnetal gap soliton 
supported by the system satisfy an anti Vakhitov-Kolokolov criterion. 

The mean-field dynamics of a cigar-shaped BEC is satisfactorily described by means of a quasi-one- 
dimensional (QID) Gross-Pitaevskii (CP) equation [22j. The CP equation for a cigar-shaped BEC in 
optical lattices is given by 

= -l^ + Vocos{2kr.s)4> + g,M'<P, (1) 

where 0(= 0(s,r)) is the order parameter which is normalized to the total number of atoms in the 
condensate and gm = 2as/a± [231 El]- The quantities Vq and stand for the strength and wave 
number of the OL respectively. In (1) we have measured length (s), time (r) in the units of a± and 
oj^. Understandably, IjJ±_ refers to frequency of the transverse harmonic confinement that produces the 
cigar-shaped BEC and a±_ = ^^/h/muj± with m, the atomic mass. 

As noted above, the mean-field atomic interaction qid in (1) can be varied spatially only by changing 
the s-wave scattering length a^. Experimentally, this is achieved by using a proper detuned laser light 
in the vicinity of the FRs. The dependence of on the spatial coordinate s is given by [HI [2S] 

«(^) = «o + ^^^, (2) 
d + f3I{s) 

where, I{s) is the intensity of the laser light and stands for the scattering length in the absence of 
the light. The quantities a and (3 are the constants which depend on the detuning parameter S of the 
laser. Clearly, such variation of a{s) will make gi£, in (1) space-dependent. For a large detuning laser 
beam having Gaussian intensity variation, gio can be written as 

9iD = 7o + liexp [-s'^/2w'^] , (3) 

where w is the width of the Gaussian profile and, 70 and 71 are the constants. Equation (1) with gm 
in (3) will provide the governing equation for a cigar-shaped BEC in OLs having spatially localized 
nonlinear interaction. 

In section 2 we present a variational formulation using a Gaussian ansatz and derive equations for the 
parameters of the trail solution with the help of Ritz optimization procedure. In section 3 we construct 
a potential model from the newtonian equation for the center of the trial solution and analyze how the 
effective potential changes with the variation of [i) the wave number of the OL and [ii) width [w) of 
the laser beam of Gaussian intensity variation. In both cases, we focus our attention on whether the 
width of the soliton exceeds the period of the lattice or not. In section 4 we consider the dynamics of the 
condensates in the effective potential. In particular, we study time evolution of the center and velocity 
of the soliton within the framework of a variational formulation and simultaneously perform a purely 
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numerical calculation to examine stability as well as to verify the result obtained from the variational 
approach. In section 5 we envisage a study on the existence and stability of fundamental gap soliton 
in a BEC with spatially localized atomic interaction. Finally, in section 6 we make some concluding 
remarks. 

2 Variational formulation 

The initial boundary value problem in (1) can be restated as the variational problem 



The Lagrangian density £ of (1) obtained from (4) is given by 

^ = - 0*00 + ^10.1' + IgisM' + v{sM'. (5) 

Here 0s = ^ and (pr = The solution of a self-focussing nonlinear Schrodinger equation is a seek 
function [26]. The shape of a Gaussian function is almost identical with a sech function. Therefore, 
in our variational formulation we adopt a Gaussian function as a trial solution [27]| of (1). To include 
the effects OL and spatial inhomogeneous nonlinearity we make parameters of the trial solution time- 
dependent and thus introduce 

0(s, r) = A(r)eM(^)'(^-^°M)'/2] ^l^so(s^so(r)))+^Hr)] ^ (6) 

as the variational ansatz for the problem. In (6), A{t), Sq{t) and $(t) represent amplitude, center 
of mass and phase of the solitonic solution and also act as variational parameters in our formulation. 
Understandably, A{t)~^ and sq{= ds/dr) stand for the width and velocity of the soliton. 

Now we shall make use of the Ritz optimization procedure [271 EH] to obtain equations for the 
variational parameters. In this optimization procedure, the first variation of the variational functional is 
made to vanish within a set of suitably chosen trial functions. Understandably, the variational functional 
for the problem considered here is the so-called effective Lagrangian density. Inserting (6) in (5) and 
then integrating the resulting equation over s from — oo to +oo we get 

(£) = ^A{r) (^8so(r)so(r) +4$(r) - 2so(r)' + 46"^^ cos (2fciSo(r)) 

+A(r)2(l + 70 v^) + ^ ' e-^+^-'M^)-' (7) 

^l/{2w^) + 2A{rY j 

as the effective Lagrangian density. From the vanishing conditions of the variational derivatives of (£) 
with respect to A{t) and Sq, we obtain the following equations. 

^[AA(r)]=0, (8a) 
°- ' V^qCos (2A;^5o) e'^M^ + — ^ ' e Wa(.)2 | (g^) 
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and 



8so(r)4so(r) + 4<i>(r) - 2sg + AVq cos{2kLSo)e 




(l + V270 



87271^3 
1 + 8^2 



1 + 4w;2A2 



Clearly, these equations show how the parameters of the lattice and a BEC soliton are related with 
each other. It is interesting to note that the equations in (8) provide a natural basis to study the 
dynamics of the solitons. For example, (86) is a newtonian equation that describes the evolution of the 
soliton's center in the presence of localized nonlinearity and an OL while (8a) represents an amplitude 
equation which shows that during evolution number of atoms N in the soliton remains constant since 
normalization condition of (6) gives y/nAlr) = N. On the other hand, (8c) describes the phase evolution 
of the condensate. In view of this, we shall judiciously make use of (8) to construct an effective potential 
for the soliton and then study their dynamics for different environment of the effective potential. 



3 Potential model 

Potential model is one of the simplest analytical models that may be used to understand dynamics of 
optical as well as BEC solitons [131 HH EZl ES]. In fact, this model is employed to calculate an effective 
potential for a moving soliton by drawing an analogy with the equation of motion of a classical particle. 
The effective potential of the soliton obtained from (86) is given by 

if - *=i 7iA(r)3 -"""ho \ 

V^ff = - Vo cos i2kLSo) e ^ + , , , \, e^+^-^-^w^ . (9) 
<3 y V2/^ + 8A{t)^ J 

For homogeneous nonlinear interaction (71 = 0) Ves in (9) gives a periodic potential. The spatial 
extension of this potentail depends on k^. On the other hand, the presence of spatially inhomogeneous 
interaction (71 7^ 0) gives rise to an effective potential (second term in Ves) the effect of which is likely 
to depend on the spatial width (w) of the (Gaussian) nonlinearity. When the spatial inhomogeneous 
atomic interaction acts in a BEC loaded in an optical lattice, then such an interaction will tend to change 
both the shape and magnitude of the trapping potential arising from the optical lattice. Therefore, the 
relative values of ki and w may be expected to have considerable effects on Ves and thus on the dynamics 
of the solitons. 

With a view to get a physical feeling from Ves on the behaviour of the soliton due to interplay between 
the optical lattice and spatial nonlinearity, we consider two cases, namely, (i) w > k^ and (ii) w < k^. 
In both (z) and (ii) we shall focus on whether A > A~'^{= soliton's width) or A(= 2'K/ki) < A~^. For 
the first case, the variation of Vefr as a function of sq for different ki values corresponding to A > A^^ 
is displayed in figure 1. In this figure, the dashed, dotted, dot-dashed and solid curves show Vc^ with 
Vo = 0.2, w = 1 and 71 = —0.5 for k^ = 0.05, 0.1, 0.3 and 0.4 respectively. The variational approach 
followed by us shows that the amplitude of the soliton remains constant with time. In view of this, we 
shall work with a constant values of A, namely, 2. From figure 1 we see that each curve for the effective 
potential consists of two distinct parts. A deep potential well around sq = and a periodic potential 
with smaller amplitudes on either sides of the well. Clearly, periodicity of the effective potential increases 
as ki increases. With the increase of /cl, Ve/f begins to deform near the well due to interplay between 
the OL and spatially localized atomic interaction. Particularly, for k^ = 0.1 there appears two peaks 
(dotted line), the amplitudes of which depend on k^. As k^ increases amplitudes of these peaks decrease 
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Figure 1: Effective potential resulting from OL and localized inhomogeneous interaction with A = 2, Vq = 
0.2,71 = -0.5, -u; = 1 and for ki = 0.05(dashed), O.l(dotted), 0.3(dot-dashed) and 0.4(solid). 



(dot-dashed line). For /cl = 0.4, amplitude of these peak (sohd curve) becomes too low that a sohton 
with sufficient energy initially situated at the minima near the well may cross the peak and will begin 

to move inside the well. 

On the other hand, the variation of V^cfr as a function of Sq with A = 0.5 for different values 
corresponding to A < A~^ displayed in figure 2. This figure shows that the curves for different /cl 
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Figure 2: Effective potential resulting from OL and localized inhomogeneous interaction with A = 0.5, Vq = 
0.2,7i = —0.5, w = 1 and for different values of fci,, namely, 0.05, 0.1, 0.3 and 0.4. Clearly, all curves merge and 
form a single curve. 

coincide and appear to form a single curve with only a deep well around Sq = 0. Clearly, the effects of 
the OL is almost negligible. As a result, we have only a potential well produced due to the localized 
nonlinear interaction. Since our aim in this paper is to study combined effects of the OL and spatially 
localized atomic interaction, we shall have to work with A < A~^. Note that the potential curve in figure 
2 has a minimum at sq = 0. So, one may expect that the soliton in this case may also be stable there. 
We have verified that the optical trap is unable to confine such solitons and thus they arc unstable. 

In view of the above, for w > kL we consider only the case A > A~^. In this case, the plot of 
as a function of Sq for different values of w is shown in figure 3. Here, values of the other parameters 
are kept same as those used in figure 1 but /cl = 0.4. In this figure, the dotted, dashed and solid curves 
give variations of Vcs for w = 0.5, 0.6 and 1 respectively. Clearly, all the curves represent periodic 
potentials with local minima at sq = 0. These local minima are separated by potential barriers from 
the minima of lattice potential. Interestingly, the depth of the intermediate (the central minimum) well 
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Figure 3: Effective potential resulting from OL and localized inhomogeneous interaction with A = 2, Vq = 
0.2,71 = -0.5, kL = 0.4 and for w = 0.5(dotted), 0.6 (dashed) and 1 (solid). 



introduced by the inhomogeneous interaction increases as w increases. But as w increases the heights 
[h) of these barriers decrease. At w = 1, the value of h becomes so small that the soliton with sufficient 
initial energy can cross the barrier and enter into the well. Therefore, it is clear that the crossover from 
a lattice minimum to the interaction energy minimum can also be taken place with the variation of 
laser-beam width. 

The potential curves we have portrayed above are the effective potentials for the soliton in a single- 
component BEC where all the atoms are in the same hyperfine state. However, it has also been possible to 
confine atoms of different hyperfine states or different species in a single trap leading to a two-component 
BEC |30]. In addition to inter-atomic interactions in a single-component BEC, the interaction between 
the atoms of different component (inter-species) in the two-component BEC has come into play. As a 
result, the effective potential for a component soliton in a coupled BEC will be affected due to this inter- 
species interaction. It is of interest to note that the effective potential in figure 1 has some similarity 
with the effective potential for the separation between of centers of the solitons in a two-component 
BEC with optical lattices [I3l El]. For example, the Ves in both cases consists of two parts, (i) a deep 
potential well and (ii) oscillatory wings on either side of the well. Interestingly, source of the periodic 
part of the effective potentials for two interacting solitons as well as for the soliton in the present 
case is same. However, origin the well in effective potentail in two cases are different. In the present 
case, this well (figure 1) appears due to presence of spatially localized inter-atomic interaction. On the 
other hand, it is shown in refs. 12 and 33 that the same well in the effective potential arises due to 
mutual interaction between the condensates. Therefore, the mutual interaction between the solitons in 
a coupled condensates may reasonably be approximated by a spatially localized interaction. Whatever 
the origin of the well in Ves, we see that the response of solitons to the variation of soliton and lattice 
parameters in both cases is almost same. For example, the solitons in both cases can be forced to move in 
localized/mutual interaction region only by varying the wave number k^. A soliton of a relatively larger 
width appears to experience an effective potential (figures 2, a potential having no oscillatory wings) 
similar to the potential that arises due to only mutual interaction between the solitons in a coupled BEC. 
Therefore, the system considered here may server to a certain extend as an useful model to get some 
realization on the behaviour of coupled solitons loaded in (linear or nonlinear) optical lattices. However, 
to confirm this remark we shall have to study dynamics of the soliton for different environments of the 
effective potential. 
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4 Dynamics of matter-wave bright solitons 



We have seen that the bright sohton in a BEC with spatially localized nonhnearity, loaded in OLs 
experiences an effective potential Ves which consists of two parts, (i) a deep well around sq — and (ii) 
oscillatory wings on either side of the well. Clearly, the well arises due to the inhomogeneous nonhnearity 
while the oscillation in V^s due to the OL. We also pointed out that the variation of wave number 
of the OL as well as width w of the Gaussian nonhnearity deforms V^g considerably only near the well 
at So — 0. Therefore, it will be an interesting curiosity to investigate how the soliton responses to such 
deformed periodic potential. To visualize this, we solve the ordinary differential equation in (86) using 
Runge-Kutta method for different values of sq and so si t — and calculate density profiles of the 
solitons. We also solve GP equation in (1) using a purely numerical routine to verify the result obtained 
from variational calculation. In particular, we shall make use of Crank-Nicholson method to solve (1) 
with the conditions 0) = sech{s — So)exp[iso{s — Sq)], 0(— 25, r) = 0(25, t) — and 70 = —1. Note 
that the numerical calculation of (1) not only serves as a supplement of the result obtained on the basis 
of a variational approach but also explicitly show whether the density profile will decay or not with time 
parameter r. 

The effective potentials (Kff) shown in figures (1) and (3) are equivalent in the sence of their shapes. 
The only exception is that in the first case the deformation in Ves arises with while in the later it is 
due to variation of w. In view of this, we consider only figure 1 (solid curve) to calculate different data 
required in calculating density profiles \tp\'^. The influence of the inhomogeneous interaction on lattice 
potentials is prominent near sq = 0. In order to visualize the effects arising due to interplay between the 
nonhnearity and lattice on the condensate, we shall focus our attention on the minimum nearest of the 
central minimum at Sq — 0. In figure 4, we display |0(s,r)p for an initial position of the soliton center 
at the local minimum of V^ff around so(0) = 3.9 having initial velocity so(0) = 0. In the same figure we 
also display variations of the soliton's center and velocity with r. In particular, figure 4(a) shows the 
plots of So and sq as a function of r while figures 4(6) and 4(c) display the evolution of density profiles 
|0(s,r)p obtain from numerical and variational calculations. Figure 4(a) clearly shows that position 
and velocity of the soliton remain unaltered with time r. On the other hand, the variation of density in 
figures 4(6) and 4(c) indicate that the sohtonic evolution is not affected by the lattice periodicity as well 
as inhomogeneous interaction. Note that analytical and numerical results are identical. The evolution 
of density profiles with unchanged amplitude for a larger scaled time r manifest that the solitons are 
stable. We have verified that this behavior of solitons is same for all other local minima of the potential. 




Figure 4: (a) Variations of position (dotted curve) and velocity (solid curve) of a soltion with r for so(0) = 3.9 
and so(0) = 0. Density profiles \(I){s,t)\'^ of solitons with r obtained from (6) variational and (c) numerical 
calculations. 



Now let us see what happens if the soliton having So(0) — 3.9 possesses very small initial velocity, 
say So(0) = 0.4, instead of zero. In this case, variations of position and velocity of the soliton is shown in 
figure 6(a) while figures 6(6) and 6(c) display evolution of density profiles with r obtained respectively 
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from variational and numerical approaches. From 5(a) we see that the position and velocity of the 




Figure 5: (a) Variations of position (dotted curve) and velocity (solid curve) of a soltion with time for so(0) = 3.9 
and so(0) = 0.1. Density profiles |^(s,r)p of solitons with r obtained from (6) variational and (c) numerical 
calculations. 



soliton rather than being constant start to vary. Clearly, the soliton first crosses the peak of the V^s 
(figure 1) and moves toward so = with increasing velocity. At sq = 0, its velocity becomes maximum. 
From this point, the soliton again starts to move up the potential slope with decreasing velocity but in 
the opposite direction. The motion continues until the soliton gets the initial velocity. This process goes 
on repeatedly. As a result, the soliton oscillates around the minimum of the potential well. The plots 
of |0(s,r)p in figures 5(a) and 5(6) show that variational calculation again agree with the numerical 
results. This figure also indicates that during evolution the soliton remains stable for a long time. 





Figure 6: (a) Variation of position (dotted curve) and velocity (solid curve) of a soltion with time for so(0) = 3.9 
and so(0) = 1. Density profiles |0(s,r)p of solitons with r obtained from (6) variational and (c) numerical 
calculations. 



It may be quite interesting to investigate the behaviour of the soliton if its initial velocity is relatively 
large, say <so(0) = 1 having same initial position. This behaviour is displayed in figure 6. Here 6(a) 
shows variation of sq and sq with r while 6(6) and 6(c) give evolution of density profiles of the solitons. 
Clearly, for such initial velocity the soliton instead of moving towards the potential well at Sq = 0, it 
moves away from sq — 3.9 with almost a constant speed. As a result, Sq only increases as r increases. 
However, while propagating the soliton shows oscillatory motion in space having periodicity same as 
that of the lattice potential. The numerical calculation (figure 6(c)) here also confirms the results in 
figure 6(6). The unattenuated propagation of the density profiles imply that solitons are stable. 

So far we have studied the dynamics of solitons with respect to a minimum near the well centred 
at So = of the effective potential. However, it remains to test the behaviour of the soliton for initial 
locations other than the minimum of effective potential. We plot in figure 7(a) variations of sq and 
So for so(0) = and so(0) = 4.5 while 7(6) and 7(c) give density profiles \ip{s,T)\'^ versus r. Figure 
7(a) clearly shows that both sq and sq vary periodically with r. However, unlike the previous case the 
periodicity of such motion depends on the displacement of the soliton from its nearest local potential 
minimum (sq = 3.9). In particular, they exhibit oscillatory motion about the nearest local potential 
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Figure 7: (a) Variation of position (dotted curve) and velocity (solid curve) of a soltion with time for sq{0) = 4.5 
and so(0) = 0. Density profiles |(/)(s,t)P of solitons with r obtained from (6) variational and (c) numerical 
calculations. 



minimum. It may be noted that the behaviour of the sohton remains same for any values of sq which 
do not correspond to potential extrema. As in the previous case, the density profiles obtained from 
variational [7(6)] and numerical [7(c)] calculations are identical and exhibit stable oscillatory motion. 



5 Gap soliton solution 

An important feature of a BEC in OLs is that it can hold a novel type of spatially localized state which 
is bright-soliton like and can exist even in a repulsive BEC [32]. This soliton resides inside a gap of 
the band-gap spectrum of linear version of the GP equation with a periodic potential and thus they 
are called gap solitons (GS). Understandably, this type of locahzed state exits due to balance between 
self-repulsion and negative effective mass provided by the OL spectrum and the chemical potential of GS 
must fall in the spectrum's band gap, where delocalized Bloch-wave states can not exist. The creation 
of gap soliton in QID BEC has also been reported in [15l |33]. Here, the BEC was pushed in a state 
with negative effective mass by means of acceleration. 

In this section, we shall consider fundamental gap solitons in a QID BEC with spatially localized 
interaction. The fundamental gap solitons [31] are indeed compact objects with quasi-Gaussian shapes, 
trapped in a single cell of the optical lattice potential. Such gap solitons in the presence of spatially 
localized interaction can be descrided by means of the GP equation (1). The stationary version of (1) 
is given by 

= + Vocos{2kLs)^ + (70 + 7iexp [-sy2w^]) l^-lV- (10) 

Here, is the chemical potential. In writing (10) we have used 0(s, r) = %p{s)e~^^'^ in (1) with 
/J*^ \i}{s)\'^ds = N. To proceed with the variational approach we first assign a Lagrangian of (10). 
One can cheek that the Lagrangian 

L = fi^' + ^ (^^^ + Vo cos{2kLs)^' + l{lo + liexp [-sy2w']) (11) 

via the Euler-Lagrange equation gives (10). In order to find a trial solution of (10), note that the 
fundamental gap soliton we considered is, in fact, have the form of a bright soliton. Therefore, the trial 
solution for the gap soliton can also be adopted as 




i;(s) = J^e-^^'\ (12) 



where variational parameters and a are respectively the amplitude and width of the soliton. Inte- 
grating the equation obtained by inserting (12) in (11) over s from —00 to +00 we get the averaged 
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Lagrangian (L) in the form 



(13) 



From the vanishing conditions of and we get 



(14) 



and 



IT a2(a2 + 4iy2)3/2 



For 7o = 7i = 0, (14) and (15) give 



and 



(15) 



(16) 



(17) 



Now, we make use of (16) and (17) to draw the band gap spectrum. The band-gap diagram for = ^ 
is displayed in figure 8. Clearly, variational method gives the well-known band-gap spectrum for the 
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Figure 8: Band - gap spectrum of the linearized GP equation in (10). In this diagram we replaced Vq by — Vq. 
The solid line between semi-infinite and first gap is obtained for even soliton while the dotted line separating 
first and second gap is drawn by considering odd solitons solution 



linearized version of (10). In order to obtain fundamental gap soliton solution we take Vq and /i values 
in the first band gap and solve (10) with the initial conditions ip{s)\s=Q = 1 and ^|s=o = 0. Here, we 
have used Vq = —5 and /i = —2. Clearly, the value of fi lies in the first gap (figure 8). The result 
obtained form the numerical integration of (1) is shown in figure 9. We are interested to fundamental 
gap solitons in a BEC with repulsive atomic interaction i.e, with positive 7 values. In view of this, 
we plot in figure 9(a) |'?/'(s)p as a function of s for 70 = 0.77 and 71 = 0.5. Here, the solid curve 
represents |V'('5)P while the dashed curve serves as a guide of the lattice potential. We have noted 
that the spatially localized atomic interaction have some effects on the effective potential for the bright 
solitons. With a view to examine effects of the localized interaction on the fundamental gap soliton we 
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Figure 9: Profiles of fundamental gap solitons for w = 1, ki = 1 witii strengtii of repulsive atomic interaction 
(a) 7o = 0.77 and 71 = 0.5 (solid curve) and (6) 70 = and 71 = 1.295 (dotted curve) 



make 70 = 0, arising background scattering length oq. In this case, the density profile obtained from 
numerical integration of (10) is portrayed in figure 9{b). The dotted curve in figure 9(6) clearly shows 
that the density profile here is also localized. Interestingly, such localized wave exists only for relatively 
higher 71 values, namely, 1.295. Therefore, a BEC in OLs with spatially localized nonlinear interaction 
may support quasi-Gaussian shaped solution. 

Up to now it is not clear whether the fundamental gap solitons discussed above are stable or not. 
In order to examine stability of such localized states we make use of Vakhitov-Kolokolov criterion [35] . 
According to this criterion localized nonlinear wave like soliton will be stable if dn/dN < 0. However, 
studies on fundamental gap soliton reveal that such solitons are linearly stable as long as dfi/dN > 
[36] . Understandably, the generalized Vakhitov-Kolokolov (VK) criterion is "inverted" for a stable 
fundamental gap soliton. This inverted VK criterion has also been explained by the approximation 
based on the averaging method [20]. Here we shall exploit variational method to justify stability of the 
profiles portrayed in figures 9(a) and 9(6) and also to verify inverted VK criterion for the gap soliton 
considered. 

From (14) and (15) we can obtain 

In figure (10) we plot dfi/dN as a function of a for w = 1 and = 1. In particular, the solid curve 
shows dn/dN in presence of both homogeneous and inhomogeneous atomic interactions (70 7^ and 
7i 7^ 0) while the dotted curve gives similar plot for 70 = but 71 7^ 0. Both the curves in this figure 
clearly show d^/dN has a positive maximum value at a = 1 (indicated by vertical line). On either sides 
of a = 1, the value of d^/dN gradually decreases. Sivan et al [10] showed that the solitons are more 
stable at the maximum value of dfi/dN. Therefore, the fundamental gap solitons considered here are 
stable if their width a ~ 1. From the figures 9(a) and 9(6), we see that the width of density profiles 
of the fundamental gap solitons obtained from direct numerical calculations have width a ^ 1 and 
thus they are stable. Understandably, the positive value of dfi/dN implies that the stable fundamental 
gap soliton in the presence of spatially localized interaction also satisfy "inverted" Vakhitov-Kolokolov 
criterion. Recently, a physical mechanism of scattering of the gap solitons in the presence of linear 
Gaussian impurity has been reported in ^7\. In particular, they observe the occurrence of repeated 
reflection, transmission, and trapping regions due to resonance of the linear impurity mode with gap 
solitons. Our study on the existence of gas soliton in the presence of nonlinear Gaussian impurity may 
motivate for further study in this direction. 
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Figure 10: Variation of dji/dN as a function of a for (i) 70 = 0.77 and 71 = 0.5 (solid curve) and (ii) 70 = 
and 71 = 1.295 (dotted curve) 

6 Conclusion 

We derived an expression of the effective potential for a soliton in BECs which were confined in an OL 
having spatially localized atomic interactions with a view to study the dynamics of the soliton. The 
localized interaction introduces a deep well at a maximum of the periodic potential due to the optical 
lattice. This well is separated from the nearest lattice minima by potential wall. The amplitude of 
which depends on lattice parameter (namely, wave number) as well as width of localized inter-atomic 
interaction. Therefore, if a soliton initially located at a minimum nearest to the well is at rest, it can 
be forced to move into the well by varying only the width of localized interaction. 

We showed that dynamical behavior of the soliton depends sensitively on its initial location in 
effective potential as well as its velocity. This is obtained by solving the appropriate GP equation by a 
purely numerical routine. We examined that, depending on the initial conditions, the density profile of 
the soltion may execute steady and oscillatory motions and can move with a constant speed without any 
attenuation. The unattcmiatcd long time evolution of the density profile implied that during evolution 
the profile remains stable. In addition, the numerical calculations not only gave validity of the method 
used but also justified the accuracy of the trial function for the model considered. 

We also showed that the effective potential resulting from the localized nonlinearity and OL for a 
soliton in a single-component BEG is almost same as that of the effective potential of two interacting 
solitons in a two-component BEG with optical lattices. Interestingly, we found that the response of the 
soliton was also similar to that of the solitons in two-component BEGs. Therefore, the effective potential 
for two interacting solitons in a BEG may be represented by a Gaussian space-dependent function. 

Finally, we noticed that the system can support fundamental gap solitons even in the absence of 
background scattering length and/or even only in the presence of spatially localized interaction. These 
solitons are stable only when they satisfy inverted Vakhitov-Kolokolov criterion. 
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